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Origin of the Low Energy Structure in Above Threshold Ionization 

A.S. Titfl and G.W.F. Drakc0 
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We present an ah initio analytic theory to account for both the very low energy structure (VLES) 
[C. Y. Wu et al., Phys. Rev. Lett. 109 , 043001 (2012); W. Quan et al, Phys. Rev. Lett. 103 , 093001 
(2009)], and the low energy structure (LES) [W. Quan et al. Phys. Rev. Lett. 103 , 093001 (2009); 
C.I. Blaga et al., Nat. Phys. 5 , 335 2009)] of above threshold ionization. The origin of both VLES 
and LES lies in a forward scattering mechanism by the Coulomb potential. We parameterize the 
S matrix in terms of a., which is the displacement of the the classical motion of an electron in 
the laser field. When a = 0, the S matrix is singular, which we attribute to be forward Coulomb 
scattering without absorption of light quanta. By devising a regularization scheme, the resulting S 
matrix is non-singular when cc = 0, and the origins of VLES and LES are revealed. We attribute 
VLES to multiple forward scattering of near-threshold electrons by the Coulomb potential, with no 
absorption of light quanta, signifying the role of the Coulomb threshold effect. We attribute LES 
to be due to the combined role of the Coulomb threshold effect and rescattering in the forward 
direction by the Coulomb potential with the absorption of light quanta. A comparison of theory 
with experiment confirms these conclusions. Further more, recently Dura et al. [Sci. Rep. 3 , 2675 
(2013)] reported the detection of slow electrons at near zero momentum, at 1.3 meV, which is much 
below the VLES, almost at threshold. Our theoretical formulation gives rise to slow electrons at 
near zero momentum and at threshold. In addition, for circularly polarized fields, it conserves the 
angular momentum in the ionization process which necessitate the disappearance of the VLES, LES 
and the slow electrons near threshold. 

PACS numbers: 42.50.Hz, 32.80.Fb, 32.80.Rm 


Above threshold ionization (ATI) is a phenomenon 
characteristic of the interaction of high-intensity lasers 
with atoms in which the atom absorbs many more pho¬ 
tons than the minimum number required to produce ion¬ 
ization. The general oscillatory features of the photoelec¬ 
tron spectrum are well understood within the strong-field 
approximation (SFA) ( Zeroth-order term in the expan¬ 
sion of the S-matrix in terms of the atomic potential) 
of Keldysh-Faisal-Reiss (KFR) l|-0, but recently, unex¬ 
pected spike-like structures 0-01 have been reported in 
the low energy region around 1 eV. This threshold struc¬ 
ture is in striking contrast to the predictions of the KFR 
theory. Furthermore, during the preparation of this pa¬ 
per, Dura et al [§] reported another surprise, the detec¬ 
tion of slow electrons at 1.3 meV, almost at threshold. 

The structures seem to be a universal feature of ATI 
in atoms and molecules. They consist of two humps: 
the first hump, called VLES, lies below 1 eV and weakly 
depends on the laser wavelength 0, 0]. The second 
hump, called LES, is characterized by a peak energy 
1 eV < Ek <3.7 eV which is typically observed using 
midinfrared lasers, and extends to higher energies (2-20 
eV) where the ATI spectra merge with the predictions 
of KFR theory 010. Beyond this, there is a plateau 
around 2U P due to direct electrons, and a higher en ergy 
plateau around 10 U p due to rescattered electrons [3-il3]. 
where U p is the ponderomotive energy. The goal of this 
paper is to present a full quantum mechanical theory 
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of the unique LES and VLES structures. They appear 
in the tunneling regime where the Kyldysh parameter 
7 = yjE b/2U p « 1, with Eb being the atomic ioniza¬ 
tion energy. In this region, one would normally expect 
the SFA to be valid. 

Various theoretical investigations have been carried 
out to understand the origin of the low energy structures. 
Numerical solution of the time dependent Schrodinger 
equation (TDSE) 0 0 provide quantitative agreement 
but little physical insight. The semiclassical model 
BH3 has been applied to explore the LES and revealed 
the essential role of the Coulomb potential in its produc¬ 
tion via forward scattering mechanism. Recently, Guo 
et al. presented an ad hoc heuristic quantum mechani¬ 
cal calculation demonstrating that the origin of the LES 
lies in the Coulomb interaction. However, the role of the 
Coulomb interaction in the production of the VLES is 
not well understood. Guo et al. attempted to account 
for both the VLES and the LES in terms of rescattering 
, but their calculations failed to display the VLES. In 
this paper, we present an ab initio analytical quantum 
mechanical formulation that simultaneously accounts for 
both the LES and the VLES. Contrary to the speculation 
of Guo al., we show that the VLES is due to Coulomb 
threshold effects via forward scattering (with no absorb- 
tion of light quanta), in accordance with Wigner thresh¬ 
old law for Coulomb attraction 0 , and the LES is due to 
Forward rescattering (with absorbtion of light quanta). 
The ab initio analytical approach does provide a trans¬ 
parent ideas of the process; thus making the understand¬ 
ing of VLES and LES complete. In addition, it does give 
rise to slow electrons near zero momentum and at thresh¬ 
old, in accordance with Wigner threshold law, consistent 
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with the findings of Dura et al., which we attribute to 
multiple Coulomb forward scattering. Furthermore, the 
Coulomb interaction is fundamental in physics, and so 
the Coulomb singularity and its regularization is of broad 
interest. 

We start with the exact expression for the time re¬ 
versed transition amplitude from a ground state (j>\ to 
a final continuum state Q (unless specified atomic 
units are used throughout) 


(5-l) fi = -i 


dt(ty f |VL <k) 


( 1 ) 


The final state ’Ff is a solution to the TDSE for an 


atomic electron interacting with a laser field 
(id t -H 0 -V L -V A )*f(r,t) = 0 


( 2 ) 


Here Va is the atomic Coulomb potential, and Vl = 
(t) ■ P + A ^l is the atom-laser interaction Hamil¬ 
tonian where P = —iS7 is the momentum operator and 
A (t) is the vector potential of the laser field. 

To first order in Va, Eq. (1) reads 


(S-l) fi «4 0) + 5 i 


:(i) 


( 3 ) 


S< 0) is the KFR direct electron term 


4 0) = - 


dt(^\t)\v L m(t)) ( 4 ) 


and is the rescattered electron term 


s£ 1} = -* 


dt f dt' (< v) (f) | VaGl +) Vh(t')cj)i (f')) 

J— oo •' —oo 

( 5 ) 

where is the Volkov wave function and G^\t, t') 

is the retarded Volkov propagator pj|. Analytical eval¬ 
uation of gives 

Sfi 1 ) = 27 t* ^2 d{E k + E B + Up — nuj) 


where Sp'* is given by Eq. ( 6 ), but with the term m = n 
excluded from the sum over m , and Sp^ corresponds to 
the m = n term in Eq. ( 6 ). 

Sp\ represents scattering by the atomic core from in¬ 
termediate continuum states with momentum q to final 
continuum states with momentum k with no exchange 
of extra photons. This represents forward scattering , 
q = k, without changing the kinetic energy of ionized 
electrons. This term does only contribute to the near 
threshold low energy direct electrons via forward scatter¬ 
ing by the Coulomb potential, which results in a spike 
in the near threshold low energy electrons. This is the 
origin of VLES. 

Sr * ^ is not only relevant to the high energy electrons 
of the ATI spectrum via backward scattering but also 
to the near threshold low energy electrons via forward 
scattering with the absorption of extra photons. Careful 
inspection of the expression for Sp^ given by Eq. ( 6 ) 
(m = n term is excluded), reveals that forward scattering 
with the exchange of one photon is significant whenever 
to = n ± 1; that is, when E q = E k ± w, and q || k. 
Scattering in the forward direction (q || k), then the 


= . Now, the 


matrix element (k I Va I q) oc ->- 

LES lies in the energy range 0.037 < E k < 0.136 (in a.u.) 
and for midinfrared laser wavelengths frequencies ui < 
.031, which makes jff- < 1 in the energy range of LES. 
Setting E k = E q ±uj and using < 1, then (k | Va | 


q) 


■2E k ^u-2E k (l^^-fg I + ...) 


' r [Ek 


» 1. The 


smaller w, the larger is its value. This forward scattering 
would be negligible had we employed a screened Coulomb 
potential. Thus we are inclined to attribute the origin of 
LES to forward scattering by the Coulomb potential with 
the exchange of photons (forward rescattering). Precise 
quantitative calculations will confirm this conclusion. 

Employing the Henneberger transformation (20] 


^(-) _ e -j/*dry L (r) $ (-) 


( 8 ) 


then Eq. (2), now reads 


x / dq <q I 0i(r)> (k | Va | q) 


Up — muj 


' Eq + Eb + Up — TO W — IT] 


f ra Qn—rn ( 6 ) 


where 77 is an infinitesimal parameter, the matrix el¬ 
ements (|) are Fourier transforms, and the functions 
fm{ q), 5 n-m(q,k) are generalized and ordinary Bessel 
functions respectively. For a Coulomb potential with ef¬ 
fective charge Z, (k | V A | q) = - yyvpppr - As 77 0+, 

the integrand in Eq. ( 6 ) is divergent for to = n at q = k. 


(see supplementary material [19j). When 77 7 ^ to, it is 
nonsingular. Thus we split Sp into a regular part Si 11 
and irregular (singular) part Sp so that 


= ^ r (1) + S[ 


(i) 


d 

i-Q t ~ H 0 - V A (r -(- a) 


= 0 


(9) 


where a = 1 dr A(t). To first order in VA(r + a), 

we have 

l$P(*)>M<PP)>+l<PP)> (10) 

where | <Ff°p)) =| % k (i)) = \ k) e~ lEkt is a plane wave 
which gives the KFR term Sp' > given in Eq. (4) and 

/ OO 

dt'Gl-\t,t')V A (r'+cx(l/)) | X M 

( 11 ) 


which gives the rescattering term S^ 1 given in Eqs. (5-6). 
(7) Here Gp(f, t') is the advanced free particle propagator. 
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Now setting a = 0 (i.e. when the electron is in the 
vicinity of the atomic core) in Eq. (11) and using Eq. 
( 8 ) then the resulting S-matrix is exactly the singularity 
inm = n term in Eq. ( 6 ); i.e., S¥ (see supplementary 
material 0 ). Therefore, the singularity in the S'-matrix 
is identified to be due to a single forward scattering of 
near threshold electrons by the Coulomb potential. If 
we denote the singular component of | by | 3?;^), 

obtained from Eq. (11) by setting a = 0, then according 
to Botero and Macek 0, | can be replaced by the 

nonsingular component - |a=o> where | 4'^ are 

the Coulomb scattering states given by 


where 7 is Euler’s constant. It is clear from Eq. (13) that 
near threshold electrons suffer the greatest single forward 
scattering by the Coulomb potential. The singularity in 
Eq. (13) at k = 0 gives rise to l/k singularity in the 
ionization rates. As k —> 0, the Coulomb scattering state 

(see Eq. (12)) 4^ _) « ^ y^J 0 (2 ^{2\Zr) cos §), J 0 ( x) 
being the Bessel function. .S'-matrix with -^= singularity, 
gives a non-singular finite ionizations rates at threshold. 
Thus, the inclusion of Coulomb potential to all orders 
with multiple forward scattering is required. 

Let us define W(r, a) = V\(v + a) — Va(i"). Then Eq. 
(9) reads 


I • ’> =1 Xk) e_l7r ° /2 r(l + a) M-a, 1 , -i(kr + k • r)) 

(12) 

with a = iXZ/k, and A is a small perturbation param¬ 
eter scaling the Coulomb potential. This component 
represents a single forward scattering by the Coulomb 
center of threshold electrons without changing their en¬ 
ergy and hence the emergence of the VLES. To evalu¬ 
ate the S'-matrix due to this component, we consider the 
Coulomb scattering state to be an approximate eigen¬ 
state of e~ l f drV c( r ) [ 22 I and recognize that the arising 
space integral is a Nordsieck-type integral 23] to obtain 


—i 


f dt{e~'S tdTV ^- |^U=o|VL0i) 


Z r , t . Z -\- ih 
= -*fc[-W2 + 7 + M^—, 


(13) 


~ H 0 - VaM - W(r, a)^j } =0 


(14) 


The scattering states 4>f ' are given by 

/ OO 

dt'G^\t,t’)W(r,(x) | ¥^) 

(15) 

where G^ "* (t, t') is the Coulomb Green function. The 
Coulomb scattering states in Eq. (15) are given by 
Eq. (12) but with a = iZ/k. Replacing the ad¬ 
vanced Coulomb propagator with the free particle one, 
and using Eqs. ( 8 ) and (1) and defining 4 l( ff\t) = 
e -if drVi(T) ,j,(^ )^ we obtain 


/ OO /* OO nt 

dt(¥^\VL<k)+i dt dt'{¥-\t) | IT(r,-a(^))G[ J +) '^i(^ , )</'i(^ , )) (16) 

-00 J — 00 J — OO 


where we utilized e ~^ T dr VUr) G (+) p drV L (r) = q(+) 
and e - l -r tdTVL ( T ' ) W(a)e' l -l' tdTVL ^ = -W(-a). The 
resulting S-matrix given in Eq. (16) is nonsingular, and 
includes the Coulomb interaction to all orders in the final 
state wave function. The first term S ^ on the left is 
the direct electron term which includes multiple forward 
Coulomb scattering without changing the photoelectron 
energy (Coulomb Threshold effect). This term gives rise 
to slow electrons at threshold and to the VLES (see Fig. 
1). The second term is the nonsingular rescattered 
electron term which gives rise to the LES (see Fig. 1). 
When a = 0, W(r, a) = V\{r + a) — Va(f) = 0 thus 
assuring the removal of the singularity arising when a = 
0. Note that for r >> a, W{ r, a) « ~^’ r , which is just 
a short range potential. 

Analytical evaluation of S^ and , given in Eq. 
(16) is nontrivial. In a previous publication we presented 
an analytical evaluation of S^ [24j . Full inclusion of 


the Coulomb attraction in the the final state upholds 
the conservation of angular momentum in the ioniza¬ 
tion process [25j, which is important for circularly po¬ 
larized lasers. For linearly polarized lasers, when consid¬ 
ering electrons ionized along polarization direction, the 
resulting is the KFR term multiplied by the factor 
e -ia/2 + a ) eM */2 + 0 ) = = if. This 

would be obtained if we replace the Coulomb scattering 
states 4'^ with | x k )- This replacement is 

carried out in Eq. (16). 

Writing (S' - l) fi « -2m J2™= no d(nuj - E k - E B - 
Up) Tfi(n), T fi » T^ 0) + T f ( 1) where T ^ 0) , and are 

the T matrices associated with S^ \ and S^ respec¬ 
tively. According to Eq. (16), the squared absolute value 
|Tfi| 2 , diverges as l/k as k —> 0. As a result the 
differential ionization rate for the absorption of n pho¬ 
tons with momentum k along the polarization direction 
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a)fi(n) ~ 27rfc|T f ^°' ) +T f [ 1 -*| 2 3 4 5 6 7 8 9 10 is non-vanishing. Unlike the 
KFR, the differential ionization rate <Dg(n) remains finite 
at Ek = 0, in accordance with the Wigner threshold law 
for Coulomb attraction 0 - In the vicinity of threshold 
it is considerably larger than that predicted by the KFR 
theory (see Fig. 1(g)), thus giving rise to slow electrons 
at near zero momentum, consistent with the findings of 
Dura el al. 

In Fig. l(a-f), we present the results of the theoreti¬ 
cal calculations for the focally averaged low energy pho¬ 
toelectron spectra, along the polarization direction, of 
Ne [(a),(d)], Kr [(b),(e)], and Xe [(c),(f)] for the laser 
parameters used in the experiment of Wu et al. [gj] (a 
Gaussian envelope is assumed for the pulse shape). We 
used hydrogen-like wave functions and scattering states. 
The theoretical calculations, Fig. l(a-c), clearly indicate 
a VLES, all below 1 eV—at 0.5 eV for Ne which slightly 
shifts with intensity, and at 0.25 eV for both Kr and Xe, 
independent of intensity. Furthermore, a LES is shown, 
located at energies greater than corresponding laser fre¬ 
quency w (> 1 eV), which becomes noticeable in the long 
wavelength spectra (A > 1 (im) and shifts with intensity. 
These results are in excellent agreement with the exper¬ 
iment of Wu et al. @. To reveal the origin of VLES and 
LES, we present in Fig. l(d-f) the yield due to T^\ and 

Tg 1 - 1 respectively for Ne at intensity 380 TW/cm 2 , Kr at 
intensity 80 TW/cm 2 , and Xe at intensity 30 TW/cm 2 . 
Fig. l(d-f), indicates that the VLES is due to Tg 0) repre¬ 
senting multiple forward scattering with no absorption of 
light quanta, and the LES to be due to representing 
forward multiple scattering with the absorption of light 
quanta; i.e., forward rescattering. Furthermore, Tg 1 ' 1 is 
negligible until the electron energy becomes equal to the 
laser frequency w (w = 1.55 eV, 0.94 eV, and 0.69 eV for 
800 nrn, 1320 nrn, and 1800 nrn respectively), and that 
the LES is noticeable at longer wavelengths and becomes 
prominent at A > 1.5 /irn. Thus we conclude that the 
VLES is located at a position below the laser frequency 
aj and the LES lies beyond oj. In Fig. 1(g), we present 
the theoretical results for Xe at different intensities and 
wavelengths corresponding to a constant ponderomotive 
energy U p ~ 18 eV representing the laser parameters 


used in the experiments of Blaga et al. Q, and Guo et 
al. 0 as well as the KFR theory prediction for Xe at 
2300 nm and intensity 36 TW/cm 2 . The results indicate 
the invariance of the LES and the VLES with respect 
to the laser parameters if U p is held constant ( VLES is 
shown to be located at 0.25 eV, which wasn’t reported 
in refs, [ij and 0)- To exclude the possibility that the 
above results are due to focal averaging, we present in 
Fig. 2 the non-focally averaged low energy spectra of Kr 
at 1320 nm and intensity 80 TW/cm 2 and Xe at 1800 nm 
and intensity 30 TW/cm 2 . It clearly indicates a VLES, 
below 0.5 eV, at 0.25 eV for both Kr and Xe and a LES 
located at energies greater than the corresponding laser 
frequency w (> 1 eV). 

In conclusion, we present an ab initio analytical theory 
to account for both the VLES and LES in ATI and the 
role of the Coulomb potential in their production. By 
regularizing the ^-matrix, their origin is revealed. We 
attribute the VLES to multiple forward scattering of the 
near-threshold electrons by the Coulomb potential, with 
no absorption of extra light quanta, signifying a Coulomb 
threshold effect, and the LES to be due to the combined 
role of Coulomb threshold effects and forward Coulomb 
rescattering. The emergence of slow electrons at thresh¬ 
old and near zero momentum, which is observed by Dura 
et al., is in accordance of Wigner threshold law [17}. The 
theoretical results under the same laser parameters as 
used in the experiments of Wu et al. (fi), Blaga et al. 0 , 
and Guo et al. |l6] confirms these conclusions. Contrary 
to Ref. 0, the formulation presented here produced a 
rescattering term which is over-all smaller than the direct 
term. Finally, embedded in the analytical formulation, 
a conservation of angular momentum in the ionization 
process. The disappearance of the low energy structures 
when circularly or elliptically polarized light is used (1,3 , 
and the emergence of the LES at higher energy when el- 
liptically light is used [8|, is a consequence of the conser¬ 
vation of momentum in the ionization process. This will 
be thoroughly investigated in a future publication. 

Research support by the Natural Sciences and Engi¬ 
neering Research Council of Canada, SHARCNET, and 
Compute/Calcu Canada is gratefully acknowledged. 
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FIG. 1: (color online) Low energy pliofea&f.lSctti Dfittlgy spec¬ 
tra, along polarization direction, of Ne at 800 nm [(a), (d)], Kr 
at 1320 nm [(b),(e)], and Xe at 1800 nm [(c),(f))]. (g) is the 
low energy photoelectron spectra of Xe at different intensities 
and wavelengths corresponding to a constant ponderomotive 
potential U p ~ 18 eV. Green curve in (g) is KFR theory pre¬ 
diction. 
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(a) Kr 1320 r 



(b) Xe 1800 nm 



FIG. 2: (color online) Non-focally averaged low energy pho- 
toelectron energy spectra, along polarization direction, of Kr 
at 1320 nm and intensity 80 TW/cm 2 (a), and Xe at 1800 nm 
and intensity 30 TW/cm 2 (b). 
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DEMONESTRATION OF THE LOGARITHMIC DIVERGENCE OF THE S MATRIX IN THE 

FORWARD DIRECTION 

In the text we state that the m = n term in Eq. (6) in the text is singular in the forward direction and gives rise 
to the singular part s[2 ■ Setting m = n in Eq. (6) and using the energy conserving delta function, then we have 


Sv 1 = 27 n &{Ek + E b + U p — ruv) 

n=n o 

x J dq (q I </>i) (k I Va | q) 


Up — nco 

Eq E k IT] 


f n 90 


(1) 


The product (q | <pi) /„(q) <7o(qi k) is analytical function in q, while the matrix element (k | Va | q) ~ |k — q| 2 is 
singular when q = k. Now writing 


we obtain 


? (i) 


(q I <Pi) fn go = J dr F(r, k) e ,q r 


5(E k + E b + U p — nu>) (U p - nco) 


( 2 ) 


x / drF(r, k) / dq 


D *q r 


k - q| 2 {E q - E k - irj) 


(3) 


The integral over q is divergent in the forward direction and the divergence is logarithmic, reminiscent of the long 
range Coulomb interaction. Since E k = k 2 / 2 and E q = q 2 / 2 and using Feynmann two denominator integral formula 


a 


CD J Q 


dx 


o [Dx + (7(1 — x)] 2 


(4) 


so that 


S 


(i) 


y 5(E k + E b + U p - nw) (U p - ruv ) 

n=no 

J drF{T,k)£dxe^-^ jdp .^_^ ]2 


0 «p r 


( 5 ) 


with 


p = q- (1 - z)k 

/3 2 = x(xk 2 + 2irf) 






2 


By choosing the z-axis along p then the integral over the angles is straightforward and the remaining integral over p 
is evaluated using the residue theorem of complex variable theory so that 


c(l) 

°ir 


X 


OO 

y: $(Ek + Eb + Up - nw) (Up - nui) 

n=n o 


J dr F(r, k) e ik r 



g[—*k-rai a ; 2 fc 2 +2ia;r7r] 

\/x 2 k 2 + 2z:E77 


c(l) 

°ir 


as given by the above equation has a logarithmic divergence in the limit r) —> 0. To see this we write 


( 6 ) 


c(l) 

°ir 


~ ^ 6(Ek + E b + U p - nw) (U p - nw) 

n=n o 

x [ dr F(r, k) e* k ' r i / _ , 2;2 - 

J I Jo x 1 k* + lixr\ 


1 g[— ik-rT+i-y/ x 2 k 2 +22xr] r] 

da" — 


+ 


dx 


o \/ x 2 k 2 + 2ixr] 


The second integral on the right hand side of Eq. (7) is easily evaluated to be equal to — \ In tS-. 
Now we set i(kr — k ■ r)x = y so that for small we have 


(7) 


1 °° 

4 1} ~ t,J 2 6 ( Ek +Eb+u p- ( u p - 


x J drF(r,k)e ikr 

IT] 'i 


pi(fcr-k-r) e 2/ _ 1 

dy- 

a y 


— In 


2/c 2 . 


( 8 ) 


The integral over y is identified as an integral representation of the exponential integral function Ei(z) [ 2 ] and so we 
have 
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where 7 ' is Euler’s constant. In the limit r) —> 0, is logarithmically divergent. Since the divergence is an additive 
not a multiplicative then it can be subtracted in a regularization scheme and can be removed. 


ISOLATION OF THE SINGULARITY 

We will show now that setting a = 0 in Eq. (11) in the text gives the singularity in m = n term in Eq. ( 6 ); i.e, in 
First we need to state that y 0 (q, k) appearing in Eq. (1) in the supplementary material is a Bessel function of 
order 0, Jq(x), such that when q = k, its argument becomes equal to 0. Thus (jo (q, k) = 1 when q = k. Taking this 
into account we can write the expression for given in Eq. (1) in the supplementary material as (without any loss 
of generality, if the value of yo(Qi k) 7 C 1 when q = k, the final result is the same) 
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The dominant contribution to comes from the second part in Eq. (10) which is very large when q ss k and singular 
when q = k. Now 
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The principal value does not contribute and the 6 function causes q = k which makes the argument of the Bessel 
function g 0 to be very close to zero, which makes (go - 1) « 0. Thus the contribution from the first term is not only 
negligible compared to the singular second term but also very negligible on its own. 

Setting a = 0 in Eq. (11) in the text (VA(r + a) = Va(i")) and using Eqs. (8) and (1) in the text, then the resulting 
S matrix can be written as 
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Now we write the free particle propagator as 

G< +) (M') = -»0(i- t') J dq | Xk (t))( Xk (t') | 

where, 77 — > 0 + is implied by the outgoing boundary conditions. Thus the resulting S matrix now reads 
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Let us denote T(t) to be 

m = (-*) f* dt 1 (e-*f v ^ dT X Jt') I 

J — OO 

the plane wave | X ) =\ q) e~ lE U is an eigenvector of Vl with Vl | X q ) = Vl(q) | X q )- Writing 

| <Mt')> = | &> e * EBt ' 

e -*fv L (T)dT 1 x j t , )} = e -.s (q ,t') | q) 

where S(q, t) = \ dr (q + lA(r)) 2 is the semiclassical action of a free electron inside laser field. Proceeding and 
carrying out integration by parts over if we obtain 
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where f n are the Fourier components. Substituting Eq. (17) into Eq. (16) we obtain 
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Substituting Eq. (18) for T into Eq. (14) we obtain for the following expression 
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the temporal integral gives 2nS(E\ i + E b + U p — nw) and therefore we finally obtain 
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( 20 ) 


which, according to Eq. (10) in the supplementary material, is exactly the singularity in to = n term in Eq. (6); 
i.e., S\2. The fundamental point here which is of great importance is that the singularity is identified to be due to 
Coulomb forward scattering in the vicinity of the atomic core (since a. = 0) and this identification led to the isolation 
and hence the removal of the singularity. 
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